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Abstract Recently, the percolation transition has been characterized on in- 
teracting networks both in presence of interdependent and antagonistic inter- 
pH ' , actions. Here we characterize the phase diagram of the percolation transition 

pH ■ in two Poisson interdependent networks with a percentage q of antagonistic 

nodes. We show that this system can present a bistability of the steady state 
solutions, and both first, and second order phase transitions. In particular, 
we observe a bistability of the solutions in some regions of the phase space 
also for a small fraction of antagonistic interactions < q < 0.4. Moreover, 
we show that a fraction q > qc = 2/3 of antagonistic interactions is nec- 
essary to strongly reduce the region in phase-space in which both networks 
fJU ' are percolating. This last result suggests that interdependent networks are 

r^ \ robust to the presence of antagonistic interactions. Our approach can be ex- 

Q ■ tended to multiple networks, and to complex boolean rules for regulating the 

O ' percolation phase transition. 



Keywords Percolation • Antagonistic interactions • Interdependent 
networks 



1 Introduction 



Percolation fP,'?,'3^ is one of the most relevant critical phenomena [mSJinUZUHJ 
(^ ; [gifTOlfni[T2 , 13 , 14 . igilSlinilllllllinilllllll] that can be defined on a complex 

fT^ ■ network. Investigating the properties of percolation on single network reveals 
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the essential role played by the topology of the network in determining the 
network robustness [Tl[5] . Recently, large attention has been paid to the study 
of the percolation transition on complex networks and surprising new phe- 
nomena have been observed. On one side, new results have shown that the 
percolation can be retarded and sharpened by the Achlioptas process [231 
[231[55J[5B] . On the other side, it has been shown that when considering in- 
teracting networks, the percolation transition can be first order pTll^l^ . 
This last result is extremely interesting because a large variety of networks 
are not isolated but are strongly interacting |27ll28l[29l[30,31.32ll33ll34l[35ll36) . 
In these systems one network function depends on the operational level of 
the other networks. Examples of investigated interacting networks go from 
infrastructure networks as the power-grid [30] and the Internet to interacting 
biological networks in physiology |36j . Nodes in interacting networks can be 
interdependent, and in this case the function or activity of a node depends 
on the function of the activity of the linked nodes in the others networks. Re- 
cent results have shown that interdependent networks are more fragile than 
single networks [30l[3T| with serious implications that these results have on 
an increasingly interconnected world. 

Nevertheless, in interacting networks wc might also observe antagonistic 
interactions. If two nodes have an antagonistic interaction, the functionality, 
or activity, of a node in a network is incompatible with the functionality, of 
the other node in the interacting network. This new possibility [37], opens 
the way to introduce in the interaction networks antagonistic interactions 
that generate a bistability of the solutions. 

In this paper we introduce a fraction of antagonistic interactions in two 
otherwise interdependent networks and we study the interplay between in- 
terdependencies and antagonistic interactions. We consider this problem in 
a simplified settings by looking at two interacting Poisson networks. For two 
Poisson networks with exclusively interdependent interactions, the steady 
state of the percolation dynamics has a large region of the phase diagram in 
which both networks are percolating. In these system, a fraction q > q^ ~ 2/3 
of antagonistic interactions is necessary in order to significantly reduce the 
region in phase-space in which both networks are percolating. This show that 
interdependent networks display a significant robustness in presence of an- 
tagonistic interactions, and that also a minority of interdependent nodes is 
enough to sustain two percolating networks. 

The paper is structured as follow; in section II we will review the theory 
of percolation in single random networks, in section III we will review the 
theory of percolation in interdependent networks, in section IV we will char- 
acterize the percolation phase diagram of two Poisson networks with purely 
antagonistic interactions, in section V we will characterize the percolation 
phase diagram in networks with a fraction q of antagonistic nodes and a 
fraction \ — q oi interdependent nodes, finally in section VI we will give the 
conclusions. 



2 Percolation on single network 

Over the past ten years great attention has been paid to the percolation 
transition on single networks. The percolating cluster in a single Poisson 
network emerges at a second order phase transition when the average degree 
of the network is (k) = 1. Nevertheless, this result can change significantly 
for networks with different degree distributions. 

In order to solve the percolation problem in a random network with degree 
distribution p^ we make use of the generating functions Go{x), Gi{x) defined 
as in the following: 

Go(a;) = ^PfcX^ (1) 

k 

We indicate by S the probability that a node is part of the percolating cluster, 
and by S" the probability that following a link we reach a node that belongs to 
the percolating cluster. Each node of the network belongs to the percolating 
cluster of the network if at least one of its links brings to a node which is 
part of the percolating cluster of the network. Expressing this observation in 
terms of S and S", we obtain the relation 

S = [I - G„{1 - S')]. (2) 

Moreover, the probability 5' can be found by solving the following recursive 
equation valid on a locally tree-like network, 

S' = [1 - Gi{l - S')]. (3) 

These equations are the well known equations for the percolation transition 
on single network [Tl[2] with given degree distribution. Equation (|3]) has a 
non trivial solution S" > which emerges continuously at a second order 
phase transition when 



dGi{x) 



dx 



{k{k 1)) ^ ^ (4) 



(k) 

The percolating cluster will be present in the network as long as 

(fc(fc-l)) 



(fc) 



> 1- (5) 



Therefore for Poisson networks we have derived that the percolation condition 
Eq. ([5]) prescribes that the average connectivity of the network z = (k) = 
{k{k — \)) / {k) should be greater than one, i.e. we must have z > 1 for the 
network to be percolating. For scale-free networks, with power-law degree 
distribution p(k) oc k~^ , the percolation condition Eq. ^ implies that the 
network, as long as the power-law exponent 7 < 3, is always percolating in 
the thermodynamic limit N ^ 00. Indeed in this case the second moment 
of the degree distribution is diverging with the network size, i.e. (fc^) — > cxd 



for TV — >■ oo . This is a crucial result in complex networks theory and implies 
that scale-free networks with exponent 7 < 3 are more robust than any 
other network with finite second moment of the degree distribution, i.e. with 

()fc2) < 00. 



3 Percolation on two interdependent networks 

In this section we will review the theory of percolation on two interdependent 
networks following the approach developed by Son et al. [29] . We will assume 
that the two networks are called network A and network B and that both 
networks have the same number of nodes N . In other words our interacting 
networks constitute a multiplex. In fact each node is represented in both 
networks. A node i belongs to the percolating cluster of the interdependent 
networks if the two following condition are met 

— (i) at least one of the neighbour nodes of i in network A belongs to the 
percolating cluster of the interdependent networks. 

— (ii) at least one of the neighbour nodes of i in network B belongs to the 
percolating cluster of the interdependent networks. 

If we denote by S the probability that a node belongs to the percolating 
cluster of two interdependent networks and by S' the probability that fol- 
lowing a link we reach a node in the percolating cluster of the interdependent 
networks we have 

S = [l-G^{l-S'^)][l-Gi{l-S')B] (6) 

The recursive equations for 5^ and S'^, valid on a tree-like random network 
are given by 

S'^^il- Gf (1 - S'^m G^{1 S'^)]. (7) 

The interesting new result is that now the percolation transition can be also 
first order [30l[27ll28ll29] as the following paragraphs show in simple cases. 



3.1 Two Poisson networks with equal average degree z. 

The percolation on interdependent networks was first studied in |301I27) and 
then further characterized in [55] . A relevant example of interdependent net- 
works is represented by two Poisson networks with the same average degree 
z ~ {k) A = {k)B- In the case of Poisson networks the generating functions 
are given by G^{x) = G^{x) = G^{x) = Gf (x) = e^f^^-^'. Therefore we 
have a relevant simplification of our Eqs. ([H]) — © because S* = S'^ = 5^. 
The equation for S (Eqs. (O, (O) now reads 

5=[l-e-^]^ (8) 
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Fig. 1 Plot of the function giji) for different values of average connectivity z. At 
z = Zc = 2.455 ... a new non-trivial solution of the function g{S) = indicates the 



onset of a first order phase transition. 



By defining 5(5) == S*- [1 -e^^^J^ the Eq. © is equivalent to g(S) = 0. This 
equation has always the solution S* = but as a function for z = Zc the curve 
g[S) is tangential to the x axis and another non trivial solution emerges. 
The point z = Zc can be found by imposing the condition 



g{S) = 
dg{S) 



dS 



0, 



(9) 



identifying the point when the function g{S) is tangential to the x axis. 
Solving this system of equations we get z = Zc = 2.455407 . . . and Sc = 
0.511699.... In Figure [T] we show a plot of the function g{S) for different 
values of the average connectivity of the network z below and above the first 
order phase transition z = Zc. For z < Zc the only solution to Eq. (|5]) is 
5 = 0, for z = Zc a new non trivial solution emerge with S = Sc- Therefore 
at z = Zc we observe a phase transition of the first order in the percolation 
problem. 



3.2 Two Poisson networks with difi'erent average degree za 7^ zb 

Another important example of interdependent networks is the the case in- 
vestigated in [55] of two Poisson networks with different average degrees 
{k)A — ZA and (fc)^ — zb- In the case of Poisson networks the generating 
functions are given by Go(a;) = Gi{x) = e/'^^^^\ Therefore we have a rele- 
vant simplification of our Eqs. (HI) ^ © because S = S'j^ = S'g. The equation 
for S (Eqs. Q, ©) now reads 



!f(5) = 5'-(l 



^zaS 



)(1 



r^ZsS 



) = o 



(10) 
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Fig. 2 Phase diagram of two interdependent Poisson networks with average degree 
ZA and zb respectively. In region I we have 5 = 0, in region II we have 5 > and 
the critical line indicates the points where the first-order transition occurs. 



The discontinuous phase transition can be found by imposing the following 
conditions 






dS 



0. 



(11) 



In Figure [5] we plot the phase diagram of the percolation process on these 
two interdependent networks. In this phase diagram we have a large region 
(Region II) in which both networks are percolating (S > 0) and we observe 
a first order percolation phase transition on the critical line of the phase 
diagram. 



4 Percolation on two antagonistic networks 

In a recent paper, we have introduced antagonistic interactions in the per- 
colation of two interacting networks |37j . As in the case of interdependent 
networks we consider two networks of TV nodes. We call the networks, net- 
work A and network B respectively and every node i is represented in both 
networks, i.e. the networks form a multiplex. The difference with the case of 
interdependent network is that if a node i belongs to the percolating cluster 
of on one network it cannot belong to the percolating cluster of the other 
one. A node i belongs to the percolating cluster of network A (network B) if 
the following two conditions are met: 

— (i) at least one node reached by following the links incident to node i in 
network A (network B) belongs to the percolating cluster in network A 
(network B); 




4.5 5 



Fig. 3 Phase diagram of two antagonistic Poisson networks with average degree 
ZA and zb respectively. In region I the only stable solution is the trivial solution 
Sa = Sb = 0. In region II-A we have only one stable solution Sa > 0, Sb ~ 0, 
Symmetrically in region II-B we have only one stable solution Sa = 0, Sb > 0. 
On the contrary in region III we have two stable solutions Sa > 0, Sb = and 
Sa = 0, Sb > and we observe a bistability of the percolation steady state solution. 



— (ii) none of the nodes reached by following the links incident to node i in 
network B (network A) belongs to the percolating cluster in network B 
(network A). 

If we indicate by Sa{Sb) the probability that a node in network A (network 
B) belongs to the percolating cluster in network A (network B), and if we 
indicate by 5^(5'^) the probability that following a link in network A (net- 
work B) we reach a node in the percolating cluster of network A (network 
B), we have 



SA = [l-G^il-S'A)]Giil-S's) 
SB = [l-Giil-S's)]G^il-S'A) 



(12) 



At the same time, in a random graph with local tree structure the proba- 
bilities S'a and S'g satisfy the following recursive equations 



S'A = [l-Gt{l-S'A)]Gi{l-S's 
S's = [l-Gf{l-S's)]G^il-S'A 



(13) 



4.1 Two Poisson networks 



We consider the case of two Poisson networks with average connectivity 

{k)A = ZA and (fc)^ = zs- 



Ill this case, the generating functions take the simple expression Gi{x) = 
G^{x) = e-^^(i-^' and Gf{x) = G^{x) = e-^«(i-^). Therefore, taking 
into consideration Eas. ([T^ and Eqs. p^ we have 5^ = Sa and S'g = Sb- 
Moreover Eas. ([T5)) take the following form: 

Sa = (l-e-^-^^^)e~^^^« 

5b = (l-e-"«^«)e-"-^^-^. (14) 

These equations have always the trivial solution Sa — 0, Sb — but de- 
pending on the value of the average connectivity in the two networks, za,zb, 
other non trivial solutions might emerge. In the following we characterize the 
phase diagram described by the solution to the Eqs. (J16p keeping in mind 
that in order to draw the phase diagram of the percolation problem wc should 
consider only the stable solutions of Eqs. (fT^ as we have widely discussed in 
[57] . Here we summarize the phase diagram in Figure |31 

— Region I za < ^,zb < 1- In this region there is only the solution Sa = 
0,5i3 =0 to the Eqs. ((TBI). 

— Region II- A za > 1, zb < ln(zyi)/(l — 1/ za)- In this regions there is only 
one stable solution to the percolation problem Sa > OSb = 

— Region II-B zb > 1, za < ln(2B)/(l — 1/zb)- In this regions there is only 
one stable solution to the percolation problem Sa = OS's > 

— Region III za > 1ii(2:b)/(1 — l/^s) and zb > ln(zA)/(l — 1/za)- In this 
region we observe two stable solutions of the percolation problem with 
Sa > 0, Sb = and Sa ^ 0,Sb > 0. Therefore in this region we observe 
a bistability of the percolation configurations. 

Wc observe that in this case for each steady state configurations, only one of 
the two networks can be percolating also in the region in which we observe 
a bistability of the solutions. 



5 Percolation on interdependent networks with a fraction q of 
antagonistic nodes 

In this section wc explore the percolation phase diagram when we allow for 
a combination of antagonistic and interdependent nodes. As in the previous 
case we consider two networks of N nodes. We call the networks, network A 
and network B respectively and every node i is represented in both networks. 
If we indicate by S'a(5'b) the probability that a random node in network 
A (network B) belongs to the percolating cluster in network A (network 
B), and if we indicate by 5^(5*^) the probability that following a link in 
network A (network B) we reach a node in the percolating cluster of network 
A (network B), we have 

Sa 



Q[l- 


G^(l- 


-S'a)\G^{I^S'b) + 




+ (1- 


-Q)[l- 


-G^{1-S'a)][1-G^{1- 


-S'b)] 


?[1- 


G^il 


-S'b)]G^{1-S'j,) + 




+(1- 


-<D[i- 


-G^{1-S'b)][1-G^{1- 


-S'a)] 



(15) 




Fig. 4 Phase diagram two Poisson interdependent networks with a fraction q — 0.3 
of antagonistic interactions. 



In the same time, in a random networks with local tree structure the 
probabilities S'^ and S'g satisfy the following recursive equations 



Da 



S'^ 



q[l-Gt{l-S'^)]G^{l^S's) + 
+ {l-q)[l-Gt{l-S'Ml-G^{l 
q[l-Gf{l-S's)]G^{l-S'^) + 
+{l-q)[l-Gf{l-S'sm-G^{l 



S'b)1 
S'a)]- 



(16) 



5.1 Two Poisson networks 

We will consider the case of two interacting Poisson networks with average 
connectivities za = {k)A and zb = {k)B- We have seen that for the case of 
two fully antagonistic Poisson networks the stable percolation configurations 
correspond to states in which either one of the two networks is percolating. 
Therefore with purely antagonistic interactions the system is not able to sus- 
tain the coexistence of two percolating clusters present in both networks. 
Here we want to generalize the above case to two interacting networks with 
only a fraction q of antagonistic interactions. For two Poisson networks we 
have G^{x) = G^{x) = e^-*(=^-i) and G^{x) = Gf (cc) = e^«(^-i) and there- 
fore Sa = S'ji^ and Sb = S'g. The Eqs. (fTO)) . (fT5)) can be explicitly written in 
terms of the average connectivities of the two networks za , zb as 



Sa = fA{SA, Sb) = 
= (1 -e-^-^^'^)[{2q 



l)e- 



jSb 



+ l-q] 



10 



Region 1 
Region II-A 
Region II-B 
Region III 
Region IV 
Region V-A 
Region V-B 



Sa = Sb =0 

Sa>0,Sb=0 

Sa = 0,Sb>0 

Sa>0,Sb>0 

Sa = Sb =0 and Sa>0,Sb>0 

Sa>0-Sb ^0 and Sa>0,Sb>0 

Sa = 0,Sb>0 and Sa>0,Sb >0 



Table 1 Stable phases in the different regions of the phase diagram of the perco- 
lation on two antagonistic Poisson networks with a fraction q = 0.3 of antagonistic 
nodes (Figure |4]) 
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Fig. 5 (Color online) Hysteresis loop for q = 0.3. The hysteresis loop is performed 
using the method explained in the main text. The value of the parameter e used in 
this figure is e = 10"^. In panel (a) and (b) zb = 4.0. In panel (c) and (d) zb — 2.8. 



Sb — !b{,Sa,Sb) 



(1 



-zbSe 



)[(2g-l)e- 



-^aSa 



+ 1-9] 



(17) 



The solutions to the recursive Eqs. (flT)) can be classified into three cate- 
gories: 

— (i) The trivial solution in which neither of the network is percolating 
Sa = Sb^ 0. 

— (ii) The solutions in which just one network is percolating. In this case 
we have either S'a > 0, S's = or Sa = 0,Sb > 0. From Eqs. ^ we 



11 



Region 1 


Sa = 6's = 


Region II-A 


Sa>0,Sb =0 


Region II-B 


Sa=0,Sb >0 


Region III 


Sa>0,Sb>0 



Table 2 Stable phases in the different regions of the phase diagram of the percola- 
tion on two antagonistic Poisson networks with a fraction q — 0.45 of antagonistic 
nodes (Figure [5]). 



find that the solution Sa > 0, Sb = emerges at a critical line of second 
order phase transition, characterized by the condition 



1 

ZA = - 

q 



(18) 



Similarly the solution Sb > 0,Sa ^ emerges at a second order phase 
transition when we have 

ZB = -- (19) 

q 

Therefore we observe the phases where just one network percolates, as 
long as g > 0. This is a major difference with respect to the phase diagram 
(Figure 0) of two purelly interdependent networks. The critical lines Eqs. 
HH]) and (fTO|) are indicated as dot-dashed lines in the phase diagrams of 
the percolation transition for different value of the fraction of antagonistic 
interactions q. 

(in) The solutions for which both networks are percolating. In this case 
we have Sa > 0, Sb > 0. This solution can either emerge (a) when the 
curves Sa = /a(<S'a, Sb) and Sb = fB{SA, Sb) cross at at point Sa = ^ 
or 5*5 = (b) when the curves Sa — /a('S'a, S's) and Sb = fsiSA, Sb) 
cross at a point Sa t^ and Sb t^ where the two curves arc tangent 
one another. 

For situation (a) the critical line can be determined by imposing, for 
example, Sa -^ in Eqs. (J13p . which yields 



ZB 



tp{zA,q) 



In 



^-{1-q) /{2q-l 



q ^ 



i--{l-q) /{2q-l 



(20) 



The function i}}{zA,q) for q < 0.5 is a decreasing function of za defined 
for 1/(1 — q) < Za < 1/<Z, for q > 0.5 is an increasing function of za 
defined for l/q < za < 1/(1 — q)- For q = 0.5 the function ^{za, q) is not 
defined but has limit tpizAjq) -^ 1/g = 2. 

A condition similar to Eq. (^0)1 can be found for za, zb by using Eqs. (fTT]) 
and imposing Sb -^ 0. In particular we obtain the other critical line 



ZA = ^{zB,q)- 



(21) 



For situation (b) the critical line can be determined imposing that the 
curves Sa = Ja^Sa, Sb) and Sb = /b^Sa, Sb), arc tangent to each other 
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Fig. 6 Phase diagram two Poisson interdependent networks with a fraction q 
0.45 of antagonistic interactions. 



at the point where they intercept. This condition can be written as 



dfA 
OSa 



1 



\dSB y OSbOSa 



(22) 



where Sa, Sb must satisfy the Eqs. (|17p . This is the equation that deter- 
mines the critical hne of first-order phase transition points. 
The condition for having a tricritical point is that Ea. ipil)) or Eq. (PT|) 
are satisfied together with Eq. (|22p . If we impose that both Eq. (PU)) and 
Eq. (j22p are satisfied at the same point, the average connectivities za and 
zb must satisfy the foUowing conditions 



ZB = ip{zA,q) 
ZB = (j){zA,q) 



ZA{2q-l) 



[I -Za{1- q)][2qZA{2q - 1) + 2 - 3g] 



(23) 



If we impose that both Eq. (PT|) and Eq. ^F^ are satisfied at the same 
point, the average connectivities za and zb must satisfy the following 
conditions 



ZA = ■4^{zB,q) 
ZA = (l){zB,q) 



(24) 
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Fig. 7 Phase diagram two Poisson interdependent networks with a fraction q — 0.6 
of antagonistic interactions. 



In general the systems of Eqs. ([23)) and Eqs. ([M)) have at most two solu- 
tions each. One trivial solution to Eqs. ([25]) and Eqs. (HH) is za = zb = - 
corresponding to Sa = Sb = 0. In the following we will characterize the 
solutions to Eqs. ((23|) as a function of the fraction of the antagonist in- 
teractions q. Similar results can be drawn by studying the system of Eqs. 

dMD- 



— Case q < 0.4. The system of Eqs. (P^ has two solutions, the trivial 
solution ZA ^ Zb = - and another non-trivial solution with za < -■ 

— Case q ~ 0.4. The system of Eqs. (|23)) has only one trivial solution 
with ZA ^ Zb ^ i. Therefore the non-trivial tricritical point disappear. 

— Case 0.4 < q < 0.5. The system of Eqs. (P5|) has two solutions, the 
trivial solution za = zb = - and another non-trivial solution with 

2^ > -. It turns out that this point is not physical because it is in 
the region in which the coexistence phase Sa > and Sb > cannot 
be sustained by the system. Therefore in this region we do not have 
a non-trivial tricritical point. 

— Case 0.5 < g < | . The system of Eqs. (|23p has only the trivial solution 
ZA = Zb = -■ Therefore the non-trivial tricritical point disappear. 

— Case g > |. The system of Eqs. (P5|) has two solutions, the trivial 



solutions za = zb 



and another non-trivial solution with za > 



5.2 The phase diagram as a function of q 



As a function of the number of antagonistic interactions q the phase diagram 
of the percolation problem change significantly. In the phase diagrams re- 
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Region 1 


Sa = 6's = 


Region II-A 


Sa>0,Sb =0 


Region II-B 


Sa=0,Sb >0 


Region III 


Sa>0,Sb>0 



Table 3 Stable phases in the different regions of the phase diagram of the perco- 
lation on two antagonistic Poisson networks with a fraction q = 0.6 of antagonistic 
nodes (Figure O 




Fig. 8 Phase diagram two Poisson interdependent networks with a fraction q — 0.< 
of antagonistic interactions. 



Region I 


Sa = Sb =0 


Region II-A 


Sa>0,Sb=0 


Region II-B 


Sa = 0,Sb>0 


Region III 


Sa>0,Sb = and Sa = 0,Sb >0 


Region IV 


Sa>0,Sb>0 


Region V-A 


Sa>0,Sb=0 and Sa>0,Sb >0 


Region V-B 


Sa = 0,Sb>0 and Sa>0,Sb >0 


Region VI 


Sa> Sb>0 and Sb > Sa > 



Table 4 Stable phases in the different regions of the phase diagram of the perco- 
lation on two antagonistic Poisson networks with a fraction q = 0.8 of antagonistic 
nodes (Figure [S]). 



ported in Figures 141 61 71 81 we plot as red lines the curves along which a first 
order phase transition can be observed and as black dashed lines the critical 
lines for a second order phase transition. 

- Case q < 0.4. 

In Figurc|3]we show the phase diagram for q = 0.3 which is a typical phase 
diagram in the region < q < 0.4. The stable phases in the different 
regions of the phase space are characterized in Table [1] From this table it 
is evident that in regions IV, V-A and V-B we observe a bistability of the 
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solutions. When g — > 0, region II-A, II-B, III, V-A and V-B disappear, 
reducing the phase diagram to Figure [2J 

In order to demonstrate the bistabihty of the percolation solution in re- 
gion IV and V-A, V-B of the phase diagram we solved recursively the Eqs. 
(flTl) for zs = 4.0 (or zb = 2.8) and variable values of za (see Figure [5]). 
We start from values of za = 3, and we solve recursively the Eqs. ([TT]). 
We find the solutions Sa — Sa{za ^ 3) > 0, Sb ~ Sb{za = 3) = 0. Then 
we lower slightly za and we solve again the Eqs. (|17p recursively, starting 
from the initial condition S"^ = SAizA = 3)-l-e, Sb = SBizA = 3)-l-e, and 
plot the result. (The small perturbation e > is necessary in order not 
to end up with the trivial solution Sa = 0,Sb = 0.) Using this procedure 
we show that if we first lower the value of za and then again we raise it, 
as shown in Figure [5l the solution present an hysteresis loop. This means 
that in the region IV and V-A, V-B there is a bistabihty of the solutions. 
Case OA<q< 0.5. 

In Figure [5] we show the phase diagram for q = 0.45 which is a typical 
phase diagram in the range 0.4 < q < 0.5. The stable phases in the 
different regions of the phase space are characterized in Table [2] For this 
range of parameters we do not observe a bistabihty of the solutions. 
Case 0.5 < g < |. 

In Figure [7] we show the phase diagram for q = 0.6 which is a typical 
phase diagram in the range 0.5 < g < |. The stable phases in the differ- 
ent regions of phase space are characterized in Table |31 From this table 
it is evident that in this case we do not observe bistabihty of the solu- 
tions. Moreover from the phase diagram Figure [7] it is clear that also if 
the majority of the nodes are antagonistic the interdependent nodes are 
enough to sustain a phase in which both networks are percolating at the 
same time (Region III). 
Case g > |. 

In Figure |8] we show the phase diagram for q = 0.8 which is a typical 
phase diagram in the range g > | . In Table 2] we characterize the stable 
phases in the different regions of the phase diagram. Region III, V-A, V-B 
and VI show a bistabihty of the solutions. In Figure [S] we show evidence 
that in these regions we can observe an hysteresis loop if we proceed 
by calculating Sa, and Sb recursively from Eqs. p7[) using the same 
technique used to produce Figure [51 For g > | the regions in phase space 
where we observe the coexistence of two percolating phases (Region IV, 
V-A, V-B and VI) are reduced and disappear as 5 — > 1. 



6 Conclusions 

In this paper we have investigated how much interdependcncies and incom- 
patibilities modify the stability of complex networks and change the phase 
diagram of the percolation transition. We found that interdependent net- 
works are robust against antagonistic interactions, and that we need a frac- 
tion q > qc = 2/3 of antagonistic interactions for reducing significantly the 
region in phase-space in which both networks are percolating. Nevertheless, 
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Fig. 9 (Color online) Hysteresis loop for q = 0.8. The hysteresis loop is performed 
using the method explained in the main text. The value of the parameter e used in 
this figure is e = 10^'^. In panel (a) and (b) zb = 5.7. In panel (c) and (d) zb = 4.5. 



we observe that even a small fractions of antagonistic nodes < (j < 0.4 
might induce a bistability of the percolation solutions. In the future we plan 
to extend this model to more than two networks, including a combinatorial 
complexity [38] of dependency types to cope with the challenges of an increas- 
ingly interconnected set of technological, social and economical networks. 



References 



1. R. Cohen, K. Erez, D. Ben-Avraham, S. Havhn, Phys. Rev. Lett. 85, 4626 
(2000). 

2. M. MolUoy, B. Reed, Random Structures & Algorithms, 6 161 (1995). 

3. R. Cohen, K. Erez, D. Ben-Avraham, S. Havhn, Phys. Rev. Lett. 86, 3682 
(2001). 

4. S. N. Dorogovtsev, A. Goltsev and J. F. F. Mendes, Rev. Mod. Phys. 80, 1275 
(2008); 

5. A. Barrat, M. Barthelemy, A. Vespignani Dynamical Processes on complex Net- 
works (Cambridge University Press, Cambridge, 2008). 

6. G. Bianconi, Phys. Lett. A. 303, 166-168 (2002). 

7. S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, Phys. Rev. E 66, 
016104 (2002). 

8. M. Leone, A. Vazquez, A. Vespignani and R. Zecchina, Eur. Phys. ,J. B 28, 
191 (2002). 

9. S. Bradde, F. Caccioli, L. Dall'Asta, and G. Bianconi, Phys. Rev. Lett. 104, 
218701 (2010). 



17 



10. R. Pastor-Satorras and A. Vespignani, Phys. Rev. Lett. 86, 3200 (2001). 

11. M. A. Munoz, R. Juhasz, C. Castellano and G. Odor, Phys. Rev. Lett. 105, 
128701 (2010). 

12. S. Eubank et al. Nature 429, 180 (2004). 

13. A. E. Motter, C. Zhou and J. Kurths, Phys. Rev. E 71, 016116 (2005). 

14. T. Nishikawa, A. E. Motter, Y.-C. Lai and F. C. Hoppensteadt, Phys. Rev. 
Lett. 91, 014101 (2003). 

15. M. Barahona and L. M. Pecora Synchronization in small-world systems, Phys. 
Rev. Lett. 89, 054101 (2002). 

16. A. Arenas, A. Di'az-Guilera, J. Kurths, Y. Moreno and C. Zhou, Physics Reports 
469, 93-153 (2008). 

17. Z. Toroczkai and K. E. Sassier, Nature 428, 716 (2004). 

18. P. Echenique, J. Gomez-Gardehes, and Y. Moreno, Europhys. Lett. 71, 325 
(2005). 

19. D. De Martino, L. Dall'Asta, G. Bianconi and M. Marsih, Phys. Rev E 79, 
015101 (R) (2009). 

20. G. Bianconi, Phys. Rev. E 85, 061113 (2012). 

21. A. Halu, L. Ferretti, A. Vezzani and G. Bianconi, EPL 99, 18001 (2012). 

22. G. Bianconi, J. Stat. Mech. P07021 (2012). 

23. D. Achlioptas, R. M. D'Souza, and J. Spencer, Science 323 1453 (2009). 

24. F. Radicchi and S. Fortunato, Phys. Rev. E 81, 036110 (2010). 

25. R. A. da Costa, S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, Phys. 
Rev. Lett. 105, 255701 (2010). 

26. O. Riordan, L. Warnke, Science 333 322 (2011). 

27. J. Gao, S. V. Buldyrev, S. Havlin, and H. E. Stanley, Physical Review Letters 
107 195701 (2011). 

28. R. Parshani, S. V. Buldyrev and S. Havlin, PNAS, 108 1007 (2011). 

29. S.-W. Son, G. Bizhani, C. Christensen, P. Grassberger and M .Paczuski, EPL 
97 16006 (2012). 

30. S. V. Buldyrev, R. Parshani, G. Paul, H. E. Stanley and S. Havlin, Nature 464 
1025 (2010). 

31. A. Vespignani, Nature 464, 984 (2010). 

32. G. J. Baxter, S. N. Dorogovtsev, A. V. Goltsev and J. F. F. Mendes, Physical 
Review Letters 109, 248701 (2012). 

33. S. Gomez, A. Diaz-Guilera, J. Gomez-Gardehes, C. J. Perez- Vicente, Y. 
Moreno, A. Arenas,Phys. Rev. Lett. 110 028701 (2013). 

34. E. Cozzo, A. Arenas, Y. Moreno,Phys. Rev. E 86, 036115 (2012). 

35. J. Gomez-Gardehes, I. Reinares, A. Arenas and L.M. Floria Scientific Reports 
2 620 (2012). 

36. A. Bashan, R. P. Bartsch, J. W. Kantelhardt, S. Havlin, P. Ch. Ivanov Nature 
Communications 3 702 (2012). 

37. K. Zhao and G. Bianconi, arXiv:1210.7498l 

38. A. K. Hartmann and M. Weigt, Phase transitions in combinatorial optimization 
problems: basics, algorithms ( Wiley- VCH,Weinheim, 2005). 



